From Classical To Statistical Ocean Dynamics

Traditional ocean modeling treats fields resolved on the model grid according to the classical dynamics of continua. Variability on smaller scales is included through sundry eddy viscosities, mixing coefficients and other schemes. In this paper we develop an alternative approach based on statistical dynamics. First, we recognize that we treat probabilities of flows, not the flows themselves. Modeled dependent variables are the moments (expectations) of the probabilities of possible flows. Second, we address the challenge to obtain the equations of motion for the moments of probable flows rather than the (traditional) equations for explicit flows. For linear terms and on larger resolved scales, the statistical equations agree with classical dynamics where those of traditional modeling works well. Differences arise where traditional modeling would relegate unresolved motion to eddy viscosity, etc.. Instead, changes of entropy (<-log P> over the probability distribution of possible flows) with respect to the modeled moments act as forcings upon those moments. In this way we obtain a consistent framework for specifying the terms which, traditionally, represent subgridscale effects. Although these statistical equations are close to the classical equations in many ways, important differences are also evident; here, two phenomena are described where the results differ. We consider eddies interacting with bottom topography. It is seen that traditional eddy viscosity and/or topographic drag, which would reduce large scale flows toward rest, are wrong. The second law of thermodynamics is violated; the arrow of time is running backwards! From statistical dynamics, approximate corrections are obtained, yielding a practical improvement to the fidelity of ocean models. Another phenomenon occurs at much smaller scales in the turbulent mixing of heat and salt. Even when both heat and salt are stably stratifying, their rates of turbulent transfer should differ. This suggests a further model improvement.     

The Characteristics and Genesis of the Massive Nitrate Deposits in the Turpan-Hami Basin of Xinjiang,China

正Massive nitrate deposits have been discovered in the Turpan-Hami basin in northwestern China.Previously,large ore grade nitrate minerals were thought only to exist in the Atacama Desert in northern Chile.Estimates of the     

An example of a highly bioturbated,storm‐influenced shoreface deposit: Upper Jurassic Ula Formation,Norwegian North Sea

Integrated ichnological and sedimentological analyses of core samples from the Upper Jurassic Ula Formation in the Norwegian Central Graben were undertaken to quantify the influence of storm waves on sedimentation. Two main facies associations (offshore and shoreface) that form a progradational coarsening upward succession are recognizable within the cores. The offshore deposits are characterized by massive to finely laminated mudstones and fine‐grained sandstones, within a moderately to highly bioturbated complex. The trace fossil assemblage is dominated by deposit‐feeding structures (for example, Planolites, Phycosiphon and Rosselia) and constitutes an expression of the proximal Zoophycos to distal Cruziana ichnofacies. The absence of grazing behaviours and dominance of deposit‐feeding ichnofossils is a reflection of the increased wave energies present (i.e. storm‐generated currents) within an offshore setting. The shoreface succession is represented by highly bioturbated fine‐grained to medium‐grained sandstones, with intervals of planar and trough cross‐bedding, thin pebble lags and bivalve‐rich shell layers. The ichnofossil assemblage, forming part of the Skolithos ichnofacies, is dominated by higher energy Ophiomorpha nodosa ichnofossils and lower energy Ophiomorpha irregulaire and Siphonichnus ichnofossils. The presence of sporadic wave‐generated sedimentary structures and variability in ichnofossil diversity and abundance attests to the influence of storm‐generated currents during deposition. As a whole, the Ula Formation strongly reflects the influence of storm deposits on sediment deposition; consequently, storm‐influenced shoreface most accurately describes these depositional environments.     

Investigation of an early Pleistocene marine osmium isotope record from the eastern equatorial Pacific

Osmium isotope composition (¹⁸⁷Os/¹⁸⁸Os) and concentrations of Os, Ir and Pt are reported for an early Pleistocene section from the ODP Site 849 in the eastern equatorial Pacific. Using the data obtained in this study, the contributions from detrital and extraterrestrial particulate matter to Os concentration and ¹⁸⁷Os/¹⁸⁸Os of sediment are estimated. Our calculations show that detrital contributions to sedimentary Os are too small (<2%) to significantly shift measured bulk sediment ¹⁸⁷Os/¹⁸⁸Os away from seawater values. A moderate but significant negative correlation between ¹⁸⁷Os/¹⁸⁸Os and ³He/¹⁸⁸Os indicate that the average particulate extraterrestrial Os flux to this site is 1.21 ± 0.47 pg cm⁻² kyr⁻¹, which constitutes ?3% of total Os burial flux. The estimates of detrital and extraterrestrial Os are used to calculate the seawater ¹⁸⁷Os/¹⁸⁸Os in the early Pleistocene. The most notable features of this early Pleistocene ¹⁸⁷Os/¹⁸⁸Os record are: (1) glacial-interglacial ¹⁸⁷Os/¹⁸⁸Os differences are insignificant within errors of estimates, (2) glacial ¹⁸⁷Os/¹⁸⁸Os values are higher compared to those reported for the late Pleistocene glacials. Comparison of ¹⁸⁷Os/¹⁸⁸Os values at Site 849 to the late Pleistocene records suggests that average seawater ¹⁸⁷Os/¹⁸⁸Os change has been modest (∼5%) since the early Pleistocene. Assuming that ¹⁸⁷Os/¹⁸⁸Os difference between the glacial periods of the late and the early Pleistocene results solely from temperature dependence of weathering rates, it has been calculated that average surface temperature during the late Pleistocene glacials was 0.8 ± 0.2 °C lower than glacials in the early Pleistocene. This inference is consistent with temperature estimates based on a recent study of pCO₂ reconstruction in the Pleistocene. This observation based on limited studies of marine ¹⁸⁷Os/¹⁸⁸Os records seems to suggest that temperature played an important role in influencing chemical weathering during the Pleistocene glacials. However, more studies are needed to confirm if this temperature-weathering feedback was operational throughout the Pleistocene. A significant down core Ir-³He co-variation coupled with similar burial fluxes of Ir at Site 849 and at LL44 GPC-3 in the north Pacific point to the utility of Ir concentration as a point paleoflux tracer. However, a twofold difference in Ir burial fluxes between the eastern and the western equatorial Pacific suggests that calibration in space and time is required to use Ir concentration as a robust indicator of paleoflux through time. Significant co-variation of concentrations of Os and total alkenone during the glacials coupled with lighter δ¹³C of benthic foraminifera indicates that productivity and carbon burial played a dominant control on scavenging of Os at Site 849. In a broader context, this data set encourages future investigation of response of PGE behavior to paleoceanographic processes.     

K‐Ar age determination,whole‐rock and oxygen isotope geochemistry of the post‐collisional Bizmişen and Çaltı plutons,SW Erzincan,eastern Central Anatolia,Turkey

Post‐collisional granitoid plutons intrude obducted Neo‐Tethyan ophiolitic rocks in central and eastern Central Anatolia. The Bizmişen and Çaltı plutons and the ophiolitic rocks that they intrude are overlain by fossiliferous and flyschoidal sedimentary rocks of the early Miocene Kemah Formation. These sedimentary rocks were deposited in basins that developed at the same time as tectonic unroofing of the plutons along E–W and NW–SE trending faults in Oligo‐Miocene time. Mineral separates from the Bizmişen and Çaltı plutons yield K‐Ar ages ranging from 42 to 46 Ma, and from 40 to 49 Ma, respectively. Major, trace, and rare‐earth element geochemistry as well as mineralogical and textural evidence reveals that the Bizmişen pluton crystallized first, followed at shallower depth by the Çaltı pluton from a medium‐K calcalkaline, I‐type hybrid magma which was generated by magma mixing of coeval mafic and felsic magmas. Delta ¹⁸O values of both plutons fall in the field of I‐type granitoids, although those of the Çaltı pluton are consistently higher than those of the Bizmişen pluton. This is in agreement with field observations, petrographic and whole‐rock geochemical data, which indicate that the Bizmişen pluton represents relatively uncontaminated mantle material, whereas the Çaltı pluton has a significant crustal component. Structural data indicating the middle Eocene emplacement age and intrusion into already obducted ophiolitic rocks, suggest a post‐collisional extensional origin. However, the pure geochemical discrimination diagrams indicate an arc origin which can be inherited either from the source material or from an upper mantle material modified by an early subduction process during the evolution of the Neo‐Tethyan ocean. Copyright © 2005 John Wiley & Sons, Ltd.     

A Wind-tunnel Study of Atmospheric Boundary-Layer Flow over Vegetated Surfaces to Suppress PM10 Emission on Owens (dry) Lake

Dust storms on Owens (dry) Lake located 200 milesnortheast of Los Angeles, California, U.S.A., havecaused serious PM₁₀ emission problems. Alaboratory-based experimental study was conducted toinvestigate the efficiency of 'vegetation cover' onthe playa to suppress dust emission rates as part ofa mitigation plan.Erodible lakebed material and field vegetation (saltgrass) taken from Owens (dry) Lake were placed in awind tunnel to simulate near-surface wind flow in theatmospheric boundary layer. Vertical wind-speedprofiles and vertical PM₁₀-concentration profileswere measured over six different levels of surfacevegetation cover. In order to understand the mechanismof aeolian particle transport and dust injection intothe atmosphere, the roughness parameter, z₀, thefriction velocity, u_*, and the surface dragcoefficient, C_D, were carefully determined.An empirical equation was developed to predictPM₁₀ emission rates as a function of wind speedand level of vegetation coverage for Owens (dry) Lake.Although the current results apply only to theeffectiveness of vegetation cover to suppressPM₁₀ emissions at Owens playa, the procedure canbe employed for any arid area or region that mightutilize added vegetation as a mitigation measure.     

Reconstructing the ancestral (Plio-Pleistocene) Rio Grande in its active tectonic setting, southern Rio Grande rift, New Mexico, USA

Deposits of the ancestral Rio Grande (aRG) belonging to the Camp Rice Formation are preserved and exposed in the uplifted southern portion of the Robledo Mountains horst of the southern Rio Grande rift. The sediments are dated palaeomagnetically to the Gauss chron (upper Pliocene). The lower part of the succession lies in a newly discovered palaeocanyon cut into underlying Eocene rocks whose margins are progressively onlapped by the upper part. Detailed sedimentological studies reveal the presence of numerous river channel and floodplain lithofacies, indicative of varied deposition in channel bar complexes of low‐sinuosity, pebbly sandbed channels that traversed generally dryland floodplains and shifted in and out of the study area five times over the 1 Myr or so recorded by the succession. Notable discoveries in the deposits are: (1) complexes of initial avulsion breakout channels at the base of major sandstone storeys; (2) common low‐angle bedsets ascribed to deposition over low‐angle dunes in active channels; (3) palaeocanyon floodplain environments with evidence of fluctuating near‐surface water tables. Sand‐body architecture is generally multistorey, with palaeocurrents indicative of funnelling of initial avulsive and main fluvial discharge from the neighbouring Mesilla basin through a narrow topographic gap into the palaeocanyon and out over the study area. An avulsion node was evidently located at the stationary southern tip to the East Robledo fault during Gauss times, with aRG channels to the north flowing close to the fault and preventing fan progradation. Subsequent Matuyama growth of the fault caused (1) deposition to cease as the whole succession was uplifted in its footwall, (2) development of a thick petrocalcic horizon, and (3) fan progradation into the Mesilla basin. Parameters for the whole aRG fluvial system are estimated as: active single channels 2 m deep and 25 m wide; valley slope 0·24–0·065°; maximum mean aggradation rate 0·05 mm year^–1; major channel belt avulsion interval 200 ky; individual channel recurrence interval 100 ky; minimum bankfull mean flow velocity 1·54 m s^–1, minimum single‐channel discharge 77 m³ s^–1, bed shear stress 22·3 N m^–2; and stream power 34·3 W m^–2.     

LNAPL Recovery Endpoints: Lessons Learnt Through Modeling,Experiments, and Field Trials

It is important to estimate what light nonaqueous phase liquid (LNAPL) recovery can be practicably achieved from subsurface environments. Over the last decade, research to address this included a broad field program, laboratory measurements and experimentation, and modeling approaches. Here, we consolidate key findings from the research in the context of current literature and understanding, with a focus on a well-validated, multiphase multicomponent modeling approach to achieve estimates of reasonable endpoints for LNAPL recovery. Simple analytical models can provide approximate saturation distributions and estimates of LNAPL recoverability via transmissivity approximation, but are insufficient to predict LNAPL saturation- and composition-based recovery endpoints for various recovery technologies. This is because they cannot account for multiphase, multicomponent fate and transport and key processes such as hysteresis. Recent advances to improve estimates of the fraction of recoverable LNAPL and its transmissivity are summarized. These advances include further development and application of a well-validated model to characterize active LNAPL recovery endpoints. We present key factors that affect the determination of LNAPL recovery endpoints, and outline how recovery endpoints are affected by natural source zone depletion (NSZD—currently gaining acceptance as a LNAPL remediation option). Major factors include geo-physical characteristics of the formation, magnitude of an LNAPL release and partitioning properties of the key LNAPL constituents of concern. Based on the capabilities of the validated model, the paper also provides a basis to optimize LNAPL recovery efforts.